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ABSTRACT 

The eclipsing and double-lined spectroscopic binary V380Cyg is an extremely important 
probe of stellar evolution: its primary component is a high-mass star at the brink of leaving 
the main sequence whereas the secondary star is still in the early part of its main sequence life- 
time. We present extensive high-resolution echelle and grating spectroscopy from Ondfejov, 
Calar Alto, Victoria and La Palma. We apply spectral disentangling to unveil the individual 
spectra of the two stars and obtain new spectroscopic elements. The secondary star contributes 
only about 6% of the total light, which remains the main limitation to measuring the system's 
characteristics. We determine improved physical properties, finding masses 13.1 ± 0.3 and 
7.8 ± 0.1 Mq, radii 16.2 ± 0.3 and 4.06 ± 0.08 R©, and effective temperatures 21 750 ± 280 
and 21 600 ± 550 K, for the primary and secondary components respectively. We perform a 
detailed abundance analysis by fitting non-LTE theoretical line profiles to the disentangled 
spectrum of the evolved primary star, and reveal an elemental abundance pattern reminiscent 
of a typical nearby B star Contrary to the predictions of recent theoretical evolution models 
with rotational mixing, no trace of abundance modifications due to the CNO cycle are de- 
tected. No match can be found between the predictions of these models and the properties of 
the primary star: a mass discrepancy of 1.5 M© exists and remains unexplained. 
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1 INTRODUCTION 

The last decade has witnessed a huge improvement in modelling 
the structure and evolution of stars, particularly for higher masses. 
The inclusion of effects such as rotation and/or magnetic fields have 
caused substantial changes in the resulting predictions (see Maeder 
&. Meynet 2000; Langer et al. 2008). Some of these concern evolu- 
tionary changes in the chemical composition of stellar atmospheres. 
Due to the CNO cycle in the core of high-mass stars some elements 
are enhanced, such as helium and nitrogen, and some are depleted, 
like carbon and to a lesser extent oxygen. With deep mixing due 
to rotation, the products of core nucleosynthesis are predicted to 
be brought to the stellar surface and cause changes in atmospheric 
abundance patterns. Rotational mixing is so efficient that changes 
in the atmospheric composition should be identifiable whilst the 
star is still on the main sequence (MS). 

Anomalous abundances of helium and CNO elements have 
been noticed for many years (Leushin 1988, Lyubimkov 1998), 
and confirmed in detail by Gies & Lambert (1992), and more re- 
cently for helium abundances in B stars by Lyubimkov et al. (2004) 
and Huang & Gies (2006). However, more ambitious observational 
studies indicate that the situation is not straightforward. In a sur- 



vey of OB stars in our Galaxy and the Magellanic Clouds, multi- 
object spectroscopy for about 750 stars were secured by Evans 
et al. (2005, 2006). Elemental abundances were separately deter- 
mined for both slow (Trundle et al. 2007, Hunter et al. 2007) and 
fast (Hunter et al. 2008, 2009) rotators in the LMC. Since theoret- 
ical calculations predict the strongest effect to be the enhancement 
of the nitrogen abundance (about 0.5 dex up to the terminal age 
of the main sequence; TAMS), Hunter et al. (2008) examined the 
abundances of this element in detail. They found a rather complex 
behaviour: nitrogen enrichment is correlated with projected rota- 
tional velocity, but there are some slow rotators with a high nitro- 
gen abundance. Their findings corroborate a recent study by Morel 
et al. (2006), who analysed a sample of slowly rotating Galactic /3 
Cephei stars and found a group with nitrogen enrichment. Impor- 
tantly, Morel et al. (2006) were able to link this with enhanced mag- 
netic field strengths. It is clear that several physical phenomena and 
processes affect the surface chemical compositions of high-mass 
stars, and it is a challenge to disentangle them. 

Empirical constraints on these processes remain hard to come 
by, despite a steady improvement in observational techniques and 
capabilities (see Hilditch 2004). In this series of papers we aim to 
calibrate the abundance patterns and chemical evolution of high- 
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mass stars by analysing detached eclipsing binaries (dEBs). These 
are vital for specifying empirical constraints on the properties 
of high-mass stars, since they are the primary source of directly 
measured stellar properties (Andersen 1991). Unfortunately, accu- 
rate (2% or better) physical properties are available for only ten 
dEBs containing > 10 starfl and only four have observational 
constraints on their chemical composition. Chemical abundances 
are difficult to determine for high-mass dEBs for several reasons. 
Firstly, they tend to display only a small number of spectral lines. 
Secondly, their often high rotational velocities means the lines are 
wide and shallow. Thirdly, in dEBs the spectral lines from the two 
components interfere with each other (Tine blending'). 

In a seminal work, Simon & Sturm (1994) introduced the tech- 
nique of spectral disentangling (SPD), by which individual spectra 
of the component stars of double-lined spectroscopic binary sys- 
tems can be deduced from observations covering a range of orbital 
phases. SPD can be used to measure spectroscopic orbits which are 
not affected by line blending (see Southworth & Clausen 2007). 
The resulting disentangled spectra also have a much higher S/N 
than the original observations, so are well suited to chemical abun- 
dance analyses. As a bonus, the strong degeneracy between effec- 
tive temperature (Toff) and surface gravity (log g) is not a problem 
for dEBs because surface gravities can be measured directly and to 
high accuracy (0.01 dex or better). A detailed investigation of these 
possibilities is given by Pavlovski & Hensberge (2005; hereafter 
PH05). 

1.1 The eclipsing binary system V380 Cygni 

V380 Cyg is a very interesting probe of stellar structure and evolu- 
tion because it contains a primary (star A) which is rather evolved 
(logg = 3.1) and a less massive secondary (star B) which is not 
(log 3 = 4.1). A thorough study by Popper & Guinan (1998) 
and Guinan et al. (2000, hereafter G2000) disclosed that star A 
is the more luminous and massive component of the system, ly- 
ing near the 'blue hook' of the MS, and that star B is a ~3 mag 
fainter unevolved MS star. Due to the faintness of star B, V380 Cyg 
was for many years categorised as a single-lined spectroscopic bi- 
nary (Batten 1962 and references therein). Hill & Batten (1984) 
applied cross-correlation techniques to considerably improve the 
orbital elements. Further advances were possible only by high- 
resolution and high signal-to-noise (S/N) spectroscopy (Lyubimkov 
et al. 1996; Popper & Guinan 1998). 

The eclipses in V380Cyg are rather shallow (amplitudes of 
0.12 and 0.09 mag) and occur on an orbital period of 12. 4d. The 
radii of the stars are not determined to high accuracy despite con- 
siderable effort spent on securing a precise and complete light curve 
(G2000). These difficulties have contributed to disagreements over 
which of the two stars is hotter (Hill & Batten 1984; Lyubimkov 
et al. 1996; G2000). The best estimate of TcffA to date came 
from fitting ultraviolet and visual spectrophotometry with model 
atmosphere energy distributions, giving Tcff = 21 350 K (G2000). 
There is also a discrepancy in the derived masses of the compo- 
nents. The two most recent studies gives Ma ~ 11.1 ± 0.5 Mq 
and Mb = 6.95 ± 0.25 Mq (G2000) and Ma = 12.1 ± 0.3 Mq 
and Mb = 7.3 ± 0.3 Mq (Lyubimkov et al. 1996). 

V380 Cyg has an eccentric orbit and displays the phenomenon 
of apsidal motion with a rate of iii — 24.0 ± 1.8 degrees per 

^ An up-to-date compilation of the properties of well-studied dEBs is 
maintained at |http : / /www .astro . keele ■ ac ■ uk/~ jkt| 



lOOyr (G2000). To match the observed stellar properties to theo- 
retical evolutionary tracks these authors found that a large core- 
overshooting parameter is required, Qov = 0.6 ± 0.1. In the HR di- 
agram star A is predicted to be located near the blue point of the MS 
hook. G2000 drew a general conclusion that high-mass stars have 
large convective cores, so are more centrally condensed than pre- 
dicted by standard evolutionary theory. However, theoretical work 
by Claret (2003, 2007) does not corroborate this conclusion, find- 
ing instead that Oov = 0.41q'3 for this dEB. 

The chemical composition of V380 Cyg, in particular the he- 
lium abundances, were studied by Leushin & Topilskaya (1986) 
and Lyubimkov et al. (1996). Both groups found star A to be en- 
riched in helium and star B to be normal. A major goal of our work 
is to improve and extend these results, and thus be able to perform 
a detailed comparison with theoretical stellar models. 



2 SPECTROSCOPIC DATA 

Star B is barely visible in the optical spectrum of the binary - 
G2000 estimated the light ratio in the visual band to be La/ Lb ~ 
14.5. A successful observational programme therefore requires 
high-resolution and high-S/N spectroscopy. For this study we have 
secured spectra at four different observatories, using either diffrac- 
tion grating or echelle spectrographs. Observing logs are given in 
the Appendix (available in the electronic version of this work). 



2.1 Ondfejov spectra 

59 spectra of V380 Cyg were obtained at the Astronomical Institute 
of the Academy of Sciences of the Czech Republic in Ondfejov. 
The observations were made in the years 2004 to 2007 with the 2 m 
telescope and the coude spectrograph. The spectral interval centred 
on the Hq line covers about 500 A at a dispersion of 17 A/mm, 
whilst the spectral intervals centred on the H/9 and H7 lines cover 
around 250 A, at a dispersion 8.5 A/mm. 



2.2 Calar Alto spectra 

We obtained 43 spectra of V380 Cyg in the course of two observing 
runs (May and August 2008) at the Centro Astronomico Hispano 
Aleman (CAHA) at Calar Alto, Spain. We used the 2.2 m telescope, 
FOCES echelle spectrograph (Pfeiffer et al. 1998), and a Loral #1 li 
CCD binned 2 x 2 to decrease the readout time. With a grating angle 
of 2724, prism angle of 130 and a 150 /im slit we obtained a spec- 
tral coverage of roughly 3700-9200 A in each exposure, at a resolv- 
ing power of _R ~ 40 000. Wavelength calibration was performed 
using thorium-argon exposures, and flat-fields were obtained using 
a tungsten lamp. The observing conditions were generally good but 
several exposures were affected by clouds. 



2.3 Victoria spectra 

Spectra were obtained with the 1.2 m McKellar telescope and 
its coude spectrograph at Dominion Astrophysical Observatory 
(DAO) with a dispersion of 9 A/mm. A total of 27 spectra were 
secured in the region of the Ha line (AA6150-6760), mostly in the 
years 1996-2000. In the present work we have used only the five 
spectra secured in 2006-2007, as the orbital parameters gradually 
change due to apsidal motion. 
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2.4 La Palma spectra 

15 spectra of V380Cyg were obtained in November 2006 with the 
Nordic Optical Telescope (NOT) and FIES echelle spectrograph 
(Frandsen & Lindberg 1999) at La Palma. We used the medium- 
resolution fibre which yielded a resolving power of 7? = 47 000 
and a fixed wavelength coverage of 3640-7360 A. Wavelength cal- 
ibration was performed using thorium-argon exposures, and flat- 
fields were obtained using a halogen lamp. The observing condi- 
tions were generally reasonable (for winter in La Palma) but several 
exposures suffered from thin cloud coverage or poor seeing. 

2.5 Data reduction 

The echelle spectra (from CAHA and NOT) were bias subtracted, 
flat-fielded and extracted with the IRAI0 echelle package routines. 
Normalisation and merging of the orders was performed with great 
care, using programs developed by ourselves, to ensure that these 
steps did not causes any systematic errors in the resulting spectra. 

The Ondfejov spectra were reduced by the program SPEFO 
(Horn et al. 1996, Skoda 1996). Initial reduction of the DAO spec- 
tra (bias subtraction, flatfielding, and spectrum extraction) were 
carried out in IRAF. The wavelength calibrations were based on 
thorium-argon lamp spectra, and were performed using SPEFO. 



3 METHOD 

Our analysis follows the methods introduced by Hensberge, 
Pavlovski & Verschueren (2000) and PH05 in their studies of the 
eclipsing and double-lined spectroscopic binary V578Mon, and 
elaborated in Paper I (Pavlovski & Southworth 2009). The core of 
this approach is reconstruction of individual stellar spectra from 
the observed composite spectra using SPD, which allows a detailed 
abundance analysis using the same tools as for single stars (c.f. 
Pavlovski 2004, Hensberge & Pavlovski 2007). SPD also gives the 
velocity amplitudes of the two stars which, in combination with 
modelling of V380 Cyg's light curves, yields the physical proper- 
ties of the two stars including precise and accurate surface gravity 
measurements which are vital for the abundance analyses. 

The method of SPD returns the individual spectra of the two 
components of a binary star, but their normalisation is arbitrary. 
This is because the continuum light ratio of the two stars is not 
directly measurable from composite spectra, short of ensuring in- 
dividual spectral lines do not dip below zero flux. Whilst relative 
line strengths in a spectrum are reliable, their absolute scaling must 
be recovered from elsewhere. Here we have renormalised our dis- 
entangled spectra using the stellar light ratios found in the light 
curve analysis (Sec.O, following the procedure described in detail 
by PH05, and taking into account the line-blocking inherent to the 
Fourier disentangling method. 

We have estimated chemical abundances by fitting the renor- 
malised disentangled spectra with synthetic spectra. Non-LTE line 
formation and spectrum synthesis computations were performed 
using DETAIL and SURFACE (Giddings 1981, Butler 1984) and the 
model atoms as listed in Paper I. Model atmospheres were calcu- 
lated with ATLAS9 (Kurucz 1979). Justification for this hybrid ap- 
proach can be found in Nieva & Przybilla (2007). 

^ IRAF is distributed by the National Optical Astronomy Observatory, 
which are operated by the Association of the Universities for Research in 
Astronomy, Inc., under cooperative agreement with the NSF. 



SPD can be quite sensitive to the precise phase distribution and 
continuum normalisation of the observed composite spectra. Hynes 
& Maxted (1998) have shown that this is not a problem for radial 
velocity studies, but Hensberge, Ilijic & Torres (2008) have demon- 
strated that it is important for correct reproduction of the overall 
shape of disentangled spectra. A lack of significant time-variability 
in the relative light contributions of the two star can produce spuri- 
ous patterns in separated spectra, which in practise appear as typi- 
cally low-frequency continuum variations. Hensberge et al. (2008) 
examined this in detail and found that the phenomenon arises due 
to degeneracy of the SPD equations and/or bias progression from 
the observed spectra. Careful planning of an observing run and the 
subsequent data reduction and analysis can mostly overcome these 
difficulties, in particular by observing the target during eclipses. In 
the case of V380 Cyg we were hampered by the relative faintness of 
the secondary star, but were able to secure spectra with a very good 
distribution in orbital phase. Combined with a careful data reduc- 
tion (particularly concerning continuum normalisation and echelle 
order merging), this approach has allowed us to obtain high-quality 
disentangled spectra which can be used for abundance analysis. 



4 SPECTROSCOPIC ORBITS THROUGH SPECTRAL 
DISENTANGLING 

In the method of SPD as introduced by Simon & Sturm (1994), 
one solves both for the individual spectra of the components, and 
the optimal set of orbital parameters. Here we have performed a 
SPD analysis in Fourier space as formulated by Hadrava (1995), 
which has one big advantage. The system of equations consists of 
A'obs X A'^pix coupled equations linear in a little more than 2A'pix, 
where TVobs is the number of observed spectra, and A'^pix is the 
number of datapoints per each observed spectrum. Using discrete 
Fourier transforms this large system of equations can be coupled 
into (A*'pix/2) + 1 systems of A'obs complex equations involving 
only two unknowns in the case of binary systems, and the required 
computation time is significantly smaller. In this work we have per- 
formed SPD on a large spectral range about 3000 A wide, for which 
the speed advantage of Fourier disentangling is very important. The 
disadvantage is that weights can be assigned only to whole spectra, 
and not to individual pixels. We weighted the spectra according to 
S/N. We use the FDBINAR^code described by Ilijic et al. (2004). 

SPD was performed in spectral regions centred on the promi- 
nent He I lines. The contribution of star B to the total system light 
is rather small (La/ Lb ~ 15.5 in V), and its lines are difficult to 
see. Fortunately, star B has a similar Tcs to star A (Sec.[6l(, and so 
is best seen in the He I lines. Star B's lines also emerge in other 
regions of its disentangled spectrum, but are very weak so are of 
limited use for determining orbital elements. We therefore concen- 
trated only on the stronger helium lines. After some trial calcula- 
tions, we selected six spectral regions centred on the lines at AA 
4026, 4388, 4471 (including Mgll A4481), 4712, 4912, and 6670 
A for measuring the orbital elements. We also grouped the spec- 
tra into appropriate time intervals to account for the change in uj 
due to apsidal motion (amounting to about 6° between our first and 
last spectrum). Table[T]lists the final orbital elements, which are the 
mean and rms of the values derived from all spectral regions. 

How does our solution compare to previous results? Batten 
(1962) made a comprehensive study of all available photographic 

^ |http : / / sail ■ zpf . fer . hr/fdbinary/| 
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Table 1. Parameters of the spectroscopic orbits for V380 Cyg derived in this 
work. Throughout the SPD analysis the period was kept fixed to the value 
derived by G2000. 



Orbital period (d) 


12.425719 (fixed) 


Time of periastron passage (HJD) 


2454615.18 ±0.14 


Velocity amplitude ftTx (kms~^) 


95.1 ± 0.3 


Velocity amplitude (kms~^) 


160.5 ± 1.2 


Mass ratio q 


0.592 ±0.005 


Orbital eccentricity e 


0.206 ±0.008 


Longitude of periastron a; (°) 


134.2 ± 1.1 




- 



4465 4470 4475 4480 4485 
wavelength 

Figure 1. Observed spectra around the He I A4471 and Mgll A4481 lines. 
The upper two spectra were obtained near quadrature, and weak lines of the 
secondary are visible and indicated with a 'B'. The lower two spectra are 
disentangled spectra of the components. The faintness of star B compared 
to star A is obvious. 



spectra, and measured the lines of star B for the first time. He found 
Ka = 93.4 ± 2.0kms-^ and iCs = 160.9 ± 2.9 kms-\ which 
are very close to our values. Hill & Batten (1984) studied existing 
and new photographic spectra, and found a discrepancy between 
the systemic velocities of the two stars. Their Kb. is about 7 km s^^ 
larger than Batten's, which increases the mass of star A from 12.5 to 
13.7 M0. Lyubimkov et al. (1996) collected all available RV mea- 
surements for V380Cyg, including new CCD observation. They 
found a much smaller of 155.3 km s,~^ and a slightly larger 
of 93.95 km s^^, but gave no uncertainties. The most comprehen- 
sive investigation to date is that of Popper & Guinan (1998), based 
on high-resolution and high-S/N echelle spectra. They adopted 
the average of their own and Lyubimkov's solutions, arriving at 
Ka = 95.6 ± 0.5 km s"^ and Kb = 151 ± Skms'^ Our own 



results compare well for Ka, but a discrepancy of 10 km s exists 
for Kb. 

The main source of error in RV measurement by cross- 
correlation is the selection of the template spectrum. SPD bypasses 
this problem entirely, by directly seeking the best overall solution 
for both the spectra and the orbital elements. It therefore does not 
require guidance from a template spectrum, or pass through an in- 
termediate step of obtaining RV measurements from each observed 
spectrum. An additional advantage of this approach is that it is not 
biased by line blending. Unfortunately, a dedicated study of error 
propagation in SPD is still lacking, and the best we can do is split 
our spectra into different sets and different wavelength regions to 
obtain several independent measurements of the orbital parameters 
(Table[T]l. Recently, Southworth & Clausen (2007) have performed 
an extensive analysis of orbits derived by double Gaussian fitting, 
cross-correlation, and SPD. The SPD orbits were the most inter- 
nally consistent, followed by double Gaussian fitting. The cross- 
correlation results were affected by line blending, even after at- 
tempts were made to correct for this. Southworth & Clausen (2007) 
noticed that their SPD solutions had several local minima around the 
best fit, and used a grid search to deal with this. 

We therefore ascribe the lOkms^^ discrepancy in Kb to line 
blending in the cross-correlation solutions. Our Kb was found by 
SPD, which does not suffer from line blending or template errors, 
and should be preferred. Finally, as line blending is not a problem 
for SPD, we could use all of our spectra rather than just those where 
the lines of the two stars are resolved. Disentangled profiles of the 
He I A4471 and Mg 11 A4481 lines are shown in Fig.[T] 



5 LIGHT CURVE ANALYSIS 

Driven by suspicions that the secondary star of V380Cyg might 
be hotter than the primary star (see below), we have revisited the 
light curves presented by G2000. These data were obtained using 
the Automated Photometric Telescopes at Mt. Hopkins, Arizona, 
and total 870 observations in the UBV bandpasses. 

The light curves show substantial ellipsiodal modulation and 
reflection effect, so we have used the Wilson-Devinney (WD) code 
(Wilson & Devinney 1971; Wilson 1979, 1993), which implements 
Roche geometry and a detailed treatment of reflection and other 
physical phenomena, version of 2004/02/06. We have modified WD 
to automatically converge to the best solution using either a damped 
version of the standard differential corrections procedure (DC) or 
the robust downhill simplex algorithm AMOEBA (Press et al. 1992). 
For our final solutions we used DC and iterated until all parameter 
corrections were less than half of their formal errors. 

The full set of fixed and control parameters are given in Ta- 
ble|2l for our modelling we adopted the orbital ephemeris from 
G2000, bolometric albedos and gravity brightening exponents ap- 
propriate for radiative atmospheres (Claret 1998, Claret 2001) and 
pseudosynchronous rotation for both stars. Using the detailed re- 
flection technique of Wilson (1990) did not substantially improve 
the fit and much increased the calculation time, so was not used. 

The treatment of limb darkening can be important, and there 
is evidence that theoretically predicted coefficients are imperfect 
(Southworth 2008; Southworth et al. 2009). We therefore included 
limb darkening in three different ways: using coefficients obtained 
by bilinear interpolation in the tables of Van Hamme (1993) or 
Claret (2000), or including the coefficients as fitted parameters. The 
solutions show a negligible dependence on the way limb darkening 
is accounted for, probably because the eclipses are quite shallow 
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Table 2. Summary of the fixed and control parameters for our final WD 
solutions of the G2000 light curves of V380 Cyg. For further details please 
see the WD user guide (Wilson & Van Hamme 2004). 



Parameter 


WD name 


Star A 


StarB 


Orbital period (d) 


PERIOD 


12.425719 


Reference time (HJD) 


HJDO 


2441256.544 


Mass ratio 


RM 


0.5916 


Stellar T^ff s (K) 


TAVH, TAVC 


21500 


22000 


Rotation rates 


F1,F2 


1.0 


1.0 


Stellar albedos 


ALBl, ALB2 


1.0 


1.0 


Gravity darkening 


GRl, GR2 


1.0 


1.0 


Numerical accuracy 


Nl, n2 


60 


30 


Bolometric LD coefficient 


XBOLl, XB0L2 


0.6475 


0.6854 


U LD coefficient 


Xl, Yl 


0.3582 


0.3128 


B LD coefficient 


Xl, Yl 


0.3534 


0.3026 


V LD coefficient 


Xl, Yl 


0.3043 


0.2622 



9 p ' ' ' I ' ' ' 

B 

U X 0.7 

: xl X 5.0 



8 - 




18000 20000 22000 24000 

Teff2 (assuming Teffl=21500) 

Figure 2. Quality of the solutions for the U BV light curves individually 
and in combination, for T^ftA = 21 500 K and TcffB fixed at a range of 
values. For the individual light curves the solution quality is given by the 
rms of the residuals, and for the combined solution the reduced assuming 
observational errors equal to the rms value for each light curve. The points 
at which a fit was evaluated or the best fit was found are given by small and 
large filled circles, respectively. 



compared to the observational scatter. For our final solutions we 
used the linear limb darkening law and fixed the coefficients at the 
Van Hamme (1993) values (Table|2ll. 



5.1 The effective temperatures 

The light curves of dEBs contain almost no information about the 
absolute T^gs of the component stars: only the ratio of the T^gs is 
directly derivable from normal light curves. One of our goals was to 
investigate the quality of solutions where star B was hotter than star 
A. For one set of solutions we therefore fixed Toff a = 21 500 K 
and then fixed TcffB to values between 18 000 and 25 000 K in steps 
of 250 K. The best fit was found for each TefiB. and for the three 
light curves both individually and in combination. Fig.|2]is a plot 
of the quality of the fits for these solutions, and shows a preference 
for solutions with TcffB > Tcff a- 
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Figure 3. Comparison between the best fits found with the WD code and 
the light curves presented by G2000. The residuals are shown at the base of 
the plot with offsets from zero. In each case the top (blue) data are U, the 
middle (green) are B and the lower (red) are V. 



5.2 Final light curve solutions 

We have obtained final solutions for the UBV light curves both in- 
dividually and in combination, using our damped version of DC 
to optimise the fits. WD allows the use of predictions from Ku- 
rucz (1993) model atmospheres to link the flux ratios in different 
passbands, which means the results have a dependence on theoret- 
ical models. We prefer to avoid this option (MODE=0 and IPB = 1; 
Wilson & Van Hamme 2004), particularly as the available data are 
multi-band and include the very blue U passband. We therefore fit 
for the light contribution of each star in each passband indepen- 
dently, after assuming reasonable TcffS with which the limb dark- 
ening coefficients are derived. 

We fixed Tcsa = 21 500 K and TcSb = 22 000 K, and then 
fitted for the potentials of the two stars, the orbital inclination, ec- 
centricity and periastron longitude, a phase shift, and the light con- 
tributions of the two stars in each passband. Additional light from 
a third star was not considered (see G2000). 

The three light curves were fitted together to give the final 
parameter values, and then individually to check for consistency. 
In Table[3]we report the fitted parameters and their formal errors 
(calculated by DC from the covariance matrix), and in Fig.[3]we 
show the fits to the data. The volume-equivalent fractional radii of 
the two stars were obtained using the WD LC program. 

We caution that the formal errors can be optimistic when there 
are strong correlations between parameters, so should be consid- 
ered with care. Taking the formal errors at face value indicates poor 
agreement between light curves, with reduced values of 4.2 and 
3.4 for the stellar potentials, 3.4 for the inclination and a huge 96.9 
and 19.7 for e and uj. This highlights the limitations of formal er- 
rors, as e and lo are strongly correlated in situations such as this one 
(e.g. Southworth et al. 2004, 2007). For our final light curve param- 
eters we instead adopt the values from the combined fit to the three 
light curves and the uncertainties from the scatter between the sep- 
arate solutions for these light curves. 
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Table 3. Results of the WD code modelling process of the G2000 light curves of V380Cyg. For the final column the uncertainties come from the scatter of the 
individual solutions for the U, B and V data. For the other columns the uncertainties are formal errors calculated by the WDDC code. 



Parameter 


WD name 


U 


B 


V 


Combined 


Adopted 


otar A potential 




A f^QQ J- 07^ 


A Sdl A- A ITin 


A SS7 -1- ^'^Q 




A AOQ -t- OAf\ 


Star B potential 


PCSV 


1 n 7S -1- n 

lU. / J ^ U.ZJ 


in SI -1- n M 




1 A SI -t- n 19 


1 SI -t- 17 


Orbital inclination (°) 


XINCL 


OA 7fi -1- 

OVJ. / U ^ U.JO 


en fi3 + n 17 
ou.u J in \j. 1 1 


Q1 9< _|_ ^1 


oi 01 -1-0 91 


81 ni + 97 


Orbital eccentricity 


E 


0.2207 ± 0.0042 


0.2023 ±0.0021 


0.1831 ±0.0017 


0.1979 ±0.0017 


0.198 ±0.012 


Periastron longitude (° ) 


PERRO 


136.5 ±2.7 


140.3 ± 1.7 


150.1 ±2.1 


143.2 ±1.5 


143.2 ±5.3 


Phase shift 


PSHIFT 


-0.0459 ± 0.0009 


-0.0460 ± 0.0004 


-0.0480 ± 0.0004 


-0.0470 ± 0.0004 




U -band light from star A 


HLUM 


12.008 ±0.046 






12.055 ±0.027 


12.055 


U -band light from star B 


CLUM 


0.734 ±0.045 






0.697 ± 0.026 


0.697 


B-band light from star A 


HLUM 




9.443 ±0.023 




9.447 ± 0.025 


9.447 


B-band hght from star B 


CLUM 




0.618 ±0.022 




0.611 ±0.024 


0.611 


y-band Hght from star A 


HLUM 






9.823 ±0.024 


9.805 ±0.025 


9.805 


V-band light from star B 


CLUM 






0.619 ±0.023 


0.635 ±0.025 


0.635 


Fractional radius of star A 




0.2563 


0.2660 


0.2615 


0.2609 


0.2609 ±0.0040 


Fractional radius of star B 




0.06391 


0.06538 


0.06701 


0.06534 


0.0653 ±0.0013 


t/-band rms of residuals (mmag) 


10.386 






10.524 




B-band rms of residuals (mmag) 




5.964 




6.130 




y-band rms of residuals (mmag) 






5.782 


5.635 





6 SPECTRAL ANALYSIS OF BOTH COMPONENTS 

A model atmosphere is defined to first order by its T^g and log g, 
and when fitting observational data these two parameters can be 
quite correlated. This can be a major limition for abundance analy- 
ses of single stars, but in the case of dEBs it is possible to accurately 
measure their log gs from light and velocity curve analysis. In the 
case of V380 Cyg we have measured log — 3.136 ± 0.014 and 
log = 4.112 ± 0.017 (Secllll. 

6.1 Effective temperature determination 

The Tcff s of B stars are best determined from the silicon ionisa- 
tion balance (Becker & Butler 1990). Lines of Sill and Si III are 
present in the disentangled spectrum of V380CygA - the lack 
of Si IV lines indicates that Toff a < 24 000 K. Grids of syn- 
thetic spectra were calculated for T^e = 19000-24000 K and for 
logg — 3.136. Equivalent widths (EWs) were measured in IRAF 
for the Sill (AA 4128, 4130, 5031) and Si III (AA 4552, 4567, 
5039) lines, and calibration curves produced for Si Il/Si III EW ra- 
tios. EW measurements were made for star A and a mean value 
of Tcfi A = 21 750 ± 220 K was found from six EW ratios. Ab- 
solute EWs can only be measured from disentangled spectra once 
they have been carefully renormalised, but the ratios of EWs are 
reliable. 

A second way of measuring T^gs of B stars is to study the 
shape of the hydrogen Balmer lines. The usual degeneracy between 
Tcft and logg is not a problem here, as the logfj values of the 
stars are known, but this method does require a careful renormal- 
isation of the disentangled spectra. We used light factors derived 
from the light curve analysis (Sec.|5](, whose variation with wave- 
length is small in the region covered by US, H7 and H/3. We avoid 
Hq since it is formed high in the stellar atmosphere and is not a 
good Teff indicator. The fitting was done with a routine for multi- 
parameter optimisation, which uses a genetic algorithm to minimise 
(Tamajo et al. 2009). The projected rotational velocity was fixed 
to 11 sin i = 98 kms^^ (Sec. l6.2t . Runs were performed in which 
just Tcff (run 1) or both T^b and log g (run 2) were optimised. The 
results are tabulated in Table|4]and compared to the observed spec- 
tra in Fig.|4] 



We also tested the inclusion of the light factor as a fitted pa- 
rameter (using the separated but not renormalised spectrum) with 
encouraging results (run 3 and run 4 in Table|4}. This suggests that 
it is possible to reliably estimate TeffS using disentangled spectra 
for SB2 systems which are not eclipsing (so do not have well- 
determined light contributions from their component stars through 
eclipse analysis), which is a very interesting possibility. The results 
also corroborate the light contributions found in Sec.[5] 

The two methods agree reasonably well on the value of Tcft a, 
although the difference is slightly larger than la for run 1. Possi- 
ble sources of systematic error include the placement of the con- 
tinuum for the EW measurements (particularly for the weak Si II 
A4218 and A4231 lines), blending between Sill A4231 and Oil 
A4233, and undulations in the continuum in the region of the (wide) 
Balmer lines. Systematic errors could also arise from the adopted 
atomic data (both in the statistical equilibrium calculations and in 
the spectral synthesis) and in the LTE assumptions of spatial ho- 
mogenity, plane parallel geometry, radiative and hydrostatic equi- 
librium, etc. However, a study of B-type supergiants by Trundle et 
al. (2004), using sophisticated non-LTE analyses, yielded Tcff mea- 
surements in good agreement with the calibration of Dufton et al. 
(2000), which uses assumptions similar to those listed above. De- 
spite this, and the generally good ageement we find between the 
Tcff s obtained in different ways, we have to be cautious as it is dif- 
ficult to quantify the size of the systematic errors. We note that the 
high precison of our Tos values comes partly from having accurate 
and precise log g values for the two stars. As the final we adopted 
Tcff A ~ 21 750±280 K where error accounts for the determination 
of Toff from hydrogen line profiles, also. 



6.2 Microturbulence and projected rotational velocity 

We have estimated the microturbulence velocity, Vtuvh, by requir- 
ing that the measured abundance of a chemical element does not 
depend on the strength of its lines: the slope of a plot of EW versus 
abundance should be zero. The O II lines are the most numerous 
in the spectrum of star A, allowing us to use 32 lines for mea- 
suring titurb- We find wturb = 14 ± lkms~^, which is relatively 
high but not unusual for evolved stars with log g ^ 3 (Hunter et 
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Table 4. Results of optimised genetic algorithm fitting of the hydrogen Balmer lines of star A. In the second column the parameter which was optimised is 
indicated; parameters were otherwise fixed to log 9 = 3.131 (and £\ = 1 for the renormalised spectrum). The light curve analysis gave £\ = 0.9394. For 
comparison, analysis of the silicon ionisation balance gives T^g^ = 21 750 it 220 K. In the last column we give the difference between the value obtained in 
optimised fitting to the value adopted in this work, is also given. 



Run 


Parameter 


W 


H7 


H5 


Mean 


Difference 


1 


Teft 


21 960±85 


22 270±70 


21 860±55 


22 030±210 


280 






0.0037144 


0.0036232 


0.0038175 






2 




21 850±55 


21690±31 


21 800±45 


21 780±80 


30 




logg 


3.121±0.034 


3.125±0.031 


3.123±0.030 


3.123±0.030 


0.008 






0.0025443 


0.0023554 


0.0027782 






3 




21 830±70 


21 620±80 


21 970±75 


21 764±43 


14 




I 


0.9541±0.0021 


0.9471 ±0.0034 


0.9487±0.0029 


0.9512±0.0015 


0.0118 






0.0019254 


0.0018643 


0.0020095 






4 




21 920±95 


21 890±90 


21 845±95 


21 885±54 


135 




logs 


3.122±0.005 


3.129±0.007 


3.139±0.004 


3.136±0.003 


0.005 




I 


0.9677±0.0026 


0.9501±0.0041 


0.9675±0.0031 


0.9643±0.0018 


0.0249 




X^ 


0.0015442 


0.0015135 


0.0016024 








Figure 4. Comparison between the disentangled spectra (points) best-fitting theoretical spectra (solid lines) for H(5 (left), H7 (middle) and H/3 (right). 



al. 2007, Morel et al. 2006). It is encouraging that G2000 derived 
fturb = 12 ± lkms~^ using a completely different approach. 

The difficulty of determining uturb has been discussed re- 
cently by Hunter et al. (2008). For B stars, the rich spectrum of the 
On ion is often used for this purpose. Uncertainties can arise be- 
tween different multiplets due to errors in the adopted atomic data 
or in the strength of non-LTE effects, so using a single O II multi- 
plet is best. Dufton et al. (2005) preferred the Si III A4560 triplet, 
and this approach was adopted by Hunter et al. (2008) in a large 
VLT-FLAMES survey of B stars. We have also tried this approach, 
and found almost the same value, iiturb = 13± 1 km s~^, as above. 
The use of an alternative oxygen multiplet at 4072, 4076 and 4079 
A was less successful as the A4079 line is weak and so made very 
shallow by rotational broadening. Therefore, an error in wturbof ±2 
km s^^ would be more realistic. 

The projected rotational velocities {v sin i) of the components 
of V380 Cyg were derived from the widths of several clean spectral 
features: CII A4267, Si III A4552, Si III A4567, Oil A4591, Oil 
A4596, and O II A4662; and avoiding He I and Mg II lines (Hens- 
berge et al. 2000). We compared the observed line profiles to a set 
of theoretical spectra calculated for different v sin i values. We find 
iiA sini = 98 ± 2kms~^. The small errorbar is due to the very 
good agreement between the observed and calculated line profiles, 
and also the very high S/N (about 1000) in the disentangled spec- 
trum of star A. 



Table 5. Helium abundances determined for V380CygA from the set of 
Hel lines, adopting ritm-b = 14kms~^ as found from the Oil lines. 
The third column gives the helium abundances derived by Lyubimkov et 
al. (1996). 



Line (A) 


e(He) 


e(He) 




This work 


LRR96 


4026.2 


0.075 ± 0.003 




4387.9 


0.069 ± 0.005 


0.123 


4437.6 


0.075 ± 0.005 




4471.5 


0.099 ± 0.007 


0.153 


4713.2 


0.071 ± 0.008 


0.212 


4921.9 


0.078 ± 0.007 


0.145 


5015.7 


0.089 ±0.012 


(0.286) 


5047.7 


0.081 ± 0.008 




5875.7 


0.141 ±0.013 


(^ 0.5) 


6678.1 


0.157 ±0.018 


& 0.5) 


mean 


0.094 ± 0.031 


0.158 ±0.038 



6.3 Helium abundance 

V380CygA has been classified as spectral type B2III, in agree- 
ment with the Teff and \ogg we find here. In B2 stars lines of 
neutral helium are very strong, while lines of ionised helium are 
completely absent. We therefore calculated profiles of Hel lines 
using the non-LTE codes DETAIL and SURFACE (see Paper I for de- 
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M\\ 4712 -ITKl 4714 4715 491E1 4930 4932 4924 1926 

wavelength wavelenclh 




Figure 5. Profiles of a range of He I lines. Disentangled spectra are rep- 
resented by points, and theoretical line profiles which are the best fit to the 
disentangled profiles are shown with red lines. For comparison, profiles cal- 
culated for the solar helium abundance are also plotted with blue lines. The 
derived abundances are given in Table|5] 



tails), for helium abundances e(He) = 0.05-0.25 in steps of 0.05. 
The synthetic spectra were broadened by ?;sini = 98kms~^ and 
Vtuvh ~ 14kms~^, and e(He) was determined by minimisation 
to the disentangled and renormalised spectrum. The results for ten 
He I lines are given in Table[5] The lines of He I at A4009.3 and 
A4120.8 were not used: our model atom does not include the first 
line and the second is severely blended with O II A4121.1. The re- 
sulting mean value for the helium abundance in the photosphere of 



star A is e(He) — 0.094 ± 0.031, which is slightly above the so- 
lar value [e(He)0 = 0.089; Grevesse et al. (2007)]. The helium 
abundances found from different He I lines span a range 0.07-0.16 
(from subsolar to 80% above solar abundance) and such a large 
spread needs explaining. 

Firstly, the A4471 and A4922 line calculations are more pre- 
cise as they include transitions up to the n = A level (Przybilla 
& Butler 2001). Lyubimkov, Rostopchin & Lambert (2004, here- 
after LRL04) discuss this in detail, and also show that the other 
two diffuse He I lines (A4026 and A4388) give a systematically 
lower helium abundance. LRL04 selected the A4471 and A4922 
lines as their principal helium abundance indicators for B stars. For 
V380 Cyg A the mean helium abundance (by number of atoms) is 
A''(He) = 0.072 ± 0.006 from A4026 and A4388 whereas we find 
Af(He) = 0.089 ± 0.010 (equalling the solar value) from A4471 
and A4922. This corroborates the findings by LRL04. 

Secondly, the varying sensitivity of He I lines to f turb might 
cause abundance discrepancies. In this work we found nturb from 
O II lines, which are the most numerous lines in the spectrum of star 
A. The recent literature contains some discussion about whether 
the uturb derived from helium lines should equal that found from 
metallic lines such as O II, N II and Si III (see Hunter et al. 2007 
and references therein). LRL04 found that Wturb derived from he- 
lium lines is systematically larger than that derived from metallic 
lines. Contrary to this, during our work on V453 Cyg (Paper I) we 
simultaneously fitted the helium line profiles for abundance and 
'^turb, finding a high e(He) and a «turb substantially lower than 
those from O II lines. We rejected this possibility and instead fixed 
"VtvLih to the 14km s~^ measured from the O II lines, which resulted 
in a normal helium abundance. 

The other lines hsted in Table[5](AA 4713, 5016, 5048, 5876 
and 6678) are known to have substantial dependences on utmb- 
Moreover, the last two lines are unusually broad when comparing to 
other helium lines (Fig.|5]l which indicates that other (non-thermal) 
broadening mechanisms are present. V380CygA has evolved to 
the giant stage, and one could expect differences in the atmospheric 
velocity field compared to dwarf stars. LRL04 have discussed the 
possibility that asphericity and/or mass loss high in the atmospheres 
of hot B giants could cause excess broadening for some helium 
lines. (Herrero, Puis & Villamariz 2000). Moreover, they have es- 
timated that the A5875 and A6678 lines would give the largest he- 
lium abundances, in agreement with our findings. We will discuss 
this further in Sec.[8] 

Thirdly, it is possible that the atmosphere of star A is strati- 
fied in helium. Catanzaro (2008) has argued that variations in he- 
lium abundance derived from different lines can be explained by 
helium stratification. A comparison of his results (for the B2 V star 
HD 32123) with ours (Table|5]l reveals the same trends from line 
to line. However, a more detailed investigation of this phenomenon 
will require state-of-the-art non-LTE model atmospheres. 

6.4 CNO abundance 

The surface abundances of helium and the CNO elements are im- 
portant observational results for comparison with theoretical pre- 
dictions, because models including rotational effects predict an en- 
richment of He and N and a depletion of C and O during the evo- 
lution of B stars. The optical spectra of these objects are rich in O 
and N lines, but have fewer C lines. We determined uturb and the 
oxygen abundance in Sec. |6.2| by fitting synthetic spectra to the ob- 
served disentangled spectra of star A. Using the same approach for 
C and N gives results which show relatively low scatter between 
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Table 6. Abundances derived for V380CygA. The number of lines used 
are given in the brackets. Abundances for the solar photosphere and for 
Galactic OB stars (Hunter et al. 2009) are given in column 3, and column 
4 lists solar abundances from Grevesse et al. (2007). Abundances for the 
metallic elements are are given on the scale in which Af(He) = 12. 



Element 


Abundance 


OB stars 


Solar 


Af(He)/Af(H) 


0.098±0.010 [10] 


0.098±0.014° 


0.093±0.002 


loge(C) 


8.21±0.03 [8] 


8.00±0.19 


8.39±0.05 


loge(N) 


7.52±0.10 [24] 


7.62±0.12 


7.78±0.06 


loge(O) 


8.54±0.14 [35] 


8.63±0.16 


8.66±0.05 


loge(Mg) 


7.52±0.05 [1] 


7.25±0.17 


7.53±0.09 


loge(Si) 


7.26±0.18 [13] 


7.42±0.07 


7.51 ±0.04 


loge(Al) 


6.04±0.03 [4] 


5.94±0.14'' 


6.37±0.04 



a Derived for unevolved B stars by Lyubimkov et al. (2004) 

b Derived for B stars with T^g < 30 000 K by Daflon et al. (2003) 
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Figure 6. Comparison between the observed spectrum of star A (points) and 
the best-fitting synthetic spectrum (red line) calculated for the abundances 
listed in Table[6] The spectral range is 4625^675 A, which contains heavily 
blended lines from C II, N II and O II. 



lines. The model carbon atom we used cannot reprode the strongest 
line (C 11 A4267), but the abundance estimates for the fairly strong 
doublet on the red wing of Ha (A6556 and A6559) are consistent 
with the abundances from other C lines. 

In the spectral range 4635^655 A there is a bunch of rela- 
tively strong On, CII and Nil lines, which can provide reliable 
estimates and checks on the derived abundances for these elements 
(Pavlovski & Hensberge 2005). In Fig.|6]we show that there is good 
agreement between the disentangled spectrum of star A and a the- 
oretical spectrum calculated for the abundances derived from un- 
blended lines (Table|6j. 

N is more sensitive to changes in abundances than C, for 
which only minute changes are predicted. However, a large spread 
in N abundances have been found for MS B stars, and mecha- 
nisms other than rotational mixing may be responsible (Morel et 




4478 4480 4482 4484 

wavelength 



Figure 7. The best-fitting calculated profiles (solid line) of Mg II A4481.2, 
which is blended with the much weaker Al 111 A4479.9 line, compared to 
the observed profiles (filled symbols) for star A. 



al. 2006, Hunter et al. 2009). This is highlighted by the [N/C] ratio, 
which is a robust indicator of CNO-processed material dredged up 
to the stellar surface. Its distribution in B stars is clearly bimodal 
(Morel 2008), with one subgroup around the solar value ([N/CJq ~ 
— 0.6dex), and one with significantly higher values, ([N/C] ~ 
-0.1 dex). With [N/C] A = -0.68±0.10dex, V380Cyg A is close 
to the solar value. This is also true for the [N/O] ratio ([N/0]a = 
-1.01 ± 0.17dex and [N/O]© = -0.90dex). 



6.5 Magnesium, silicon and aluminium 

Determination of the abundances of metals like Mg, Si and Al 
makes possible an estimate of the bulk metallicity of B stars. These 
elements do not participate in the CNO cycle so their abundances 
are not affected by MS evolution. In particular, Mg is an excellent 
metallicity indicator because its solar and meteoritic abundances 
are know with high accuracy. 

A disadvantage in the use of Mg for metallicity estimates in 
B stars is that there is only one suitable spectral line in the whole 
optical region: Mg II A4481.2; although this is one of the strongest 
lines in the spectra of early-type stars. For stars with moderate and 
high i; sin i an analysis of this line is complicated by blending with 
the neighbouring AlIII A 4479.9 line. In the T^s range of interest 
here, the AlIII A4479.9 line reaches its maximum strength, and 
the ratio of the Al 111 and Mg II line EWs reaches a maximum (see 
Lyubimkov et al. 2005). 

We derived the Al abundance for star A from three lines (AA 
4512, 4671, 5766), finding a very good agreement. We used this 
measurement in calculating the Mg abundance from the Al III and 
Mg II blend at 4481 A. We find log e(Mg) = 7.52 ± 0.05, which 
translates to a metallicity index of [Mg/H] = -0.02 ± 0.05. Fig.|7] 
shows the fit of the synthetic spectrum to the data for this spectral 
line blend. 

G2000 found a subsolar metallicity for star A ([M/H] = 
—0.44 ± 0.07) by fitting the system's spectral energy distribution. 
Elemental abundances in the present work are derived from non- 
LTE line profile calculations, while G2000 used theoretical spectra 
generated in the LTE approximation. This might explain the dis- 
agreement between the two results. 
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Figure 8. Comparison of the disentangled renomialised spectrum of star B 
to tlieoretical spectra for two wavelengtli intervals, covering He I A447 1 and 
Mgll A4481 and the Si III triplet at AA4550-4580. The synthetic spectra 
were calculated for T^s= 20 000 K (blue lines), 22 000 K (red) and 24 000 
K (green). Elemental abundances were the same as for star A. The low S/N 
of the disentangled spectrum is due to relative faintness of star B. 



6.6 Secondary star 

Despite its faintness and small contribution to the total flux of the 
system, our application of SPD readily revealed the spectrum of star 
B (Fig.[T](. Its light contribution, determined from the light curve 
analysis (Sec.[5](, is 6.1% in B and 5.9% in V . To renormalise its 
disentangled spectrum to a continuum level of unity we have mul- 
tiplied it by 16.4 and 16.9, respectively. With typically 55 spectra 
available for disentangling for a given spectral range, with average 
S/N ~ 150, the reconstructed spectrum of star B has S/N ~ 45. 
This is enough to check its Tea, but too low for a detailed abun- 
dance study. 

Theoretical spectra were synthesised with Tcff = 19 000- 
25 000 K and log g = 4.112. Only strong He 1 lines, Mg II A4481 
and the Si III AA4550-4580 triplet were considered in the min- 
imisation. The mean value from eight spectral lines gives Tcff b = 
21 600 ± 550 K, which is consistent within the uncertainties with 
Tcft B and the Tcff ratio inferred from the light curve analysis. 



Table 7. The absolute dimensions and related quantities determined for 
V380Cyg. Voq and Vgynch are the observed equatorial and calculated syn- 
chronous rotational velocities, respectively. We adopt Lq = 3.826 X 
lO^e W and A/toi,© = 4.75. 





Star A 


StarB 


Semimajor axis ( R0) 


62.17 ±0.32 


Mass (M0) 


13.13 ±0.24 


7.779 ±0.095 


Radius ( R©) 


16.22 ±0.26 


4.060 ± 0.084 


log 3 (cms~^) 


3.136 ±0.014 


4.112±0.017 


Effective temperature (K) 


21 750 ±280 


21 600 ±550 


log(L/LQ) 


4.73 ± 0.028 


3.51 ±0.040 


A^bol (mag) 


-7.06 ± 0.06 


-4.03 ±0.10 


Voq (kms~l) 


98 ±2 


43±4 


Ksynch (kms-1) 


66.1 ±1.1 


16.54 ±0.34 


Extinction Eg^y (mag) 


0.20 ±0.03 


Distance using Flower V-band BCs (pc) 


1015 ± 56 


Distance using Bessell V-band BCs (pc) 


1034 ± 50 


Distance using Bessell i^T-band BCs (pc) 


1004 ± 23 


Distance using Girardi V-band BCs (pc) 


1035 ± 50 


Distance using Girardi K-hand BCs (pc) 


1004 ± 23 



In Fig.l we compare the disentangled spectrum of star B 
to theoretical spectra for two wavelength intervals, covering He I 
A4471 and Mg II A4481 and the Si III triplet. The synthetic spectra 
were generated using the He, Mg and Si abundances found for star 
A. Further progress would require a substantial increase in the S/N, 
and thus extensive new observations. 



7 THE PHYSICAL PROPERTIES OF V380 CYG 

Armed with the velocity amplitudes determined in Sec.|4]and the 
results of the light curve analysis from Sec.[5]we can determine the 
physical properties of the component stars of V380Cyg (Table|7]l. 
For this we use the JKTABSDIM code (Southworth et al. 2005), 
which propagates the uncertainties on the input quantities via a per- 
turbation analysis. The resulting error budget is dominated by the 
photometric parameters, highlighting the importance of obtaining 
a new light curve for this difficult object. Even without this, the 
masses and radii of the component stars are now known to accura- 
cies of 2% or better so are good enough for detailed tests of stellar 
evolution theory. We also find that V380 Cyg has not attained either 
orbital circularisation or rotational synchronisation. 

With the Tcff s determined above, we have found the distance 
to V380Cyg via the bolometric correction (BC) method. Empiri- 
cal BCs were obtained from Flower (1996), and theoretical ones 
from Bessell et al. (1998) and Girardi et al. (2002). We adopted 
BV magnitudes from the Tycho catalogue (H0g et al. 1997) and 
JHK magnitudes from 2MASS (Skrutskie et al. 2006). There is 
good agreement between different sources of BCs (Table|7]l. Using 
the Girardi et al. (2002) BCs, which are available for a wide range 
of passbands, we find that a reddening of Eb-v ~ 0.20 ± 0.03 is 
needed to obtain consistent distances at optical and infrared wave- 
lengths. For our final value we adopt a distance of 1004 ± 23 pc, 
from A'-band BCs. This is in good agreement with but more precise 
than the value of 1000 ± 40 pc found by G2000. 
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Figure 9. Comparison between stellar models and the absolute parameters of V380 Cyg in the log Tefp — log L and log T^g — log g diagrams. Numbers in the 
plots give the masses (in Mq) for the theoretical model tracks, and the positions of the components of V380 Cyg are indicated with filled circles. Left panels: 
models calculated with core overshooting but no rotational mixing (Schaller et al. 1992). Centre panels: models calculated by Claret (1995) for metallicity 
Z = 0.02 (blue lines) and Claret & Gimenez (1995) for Z = 0.01 (red lines), with core overshooting but no rotational mixing. Right panels: models which 
include rotational mixing (Ekstrom et al. 2008), with fi/ficrit = 0.10 (green lines), fi/ficrit = 0.50 (blue fines). Dotted red and black lines denote model 
predictions for the masses of V380 Cyg (13.1 and 7.8 Mq) for the same two rotational velocities. 



8 PROBING THE MODELS 

The two components of V380Cyg have very different masses and 
radii, and thus evolutionary stages, so are a very good test of theo- 
retical stellar models. The accuracies of the mass and radius mea- 
surements is generally 2%, limited by the quality of the available 
light curves and the faintness of the secondary star in our spectra. 
A comparison with models is helped by our abundance analysis. 
Of particular interest is the effect of rotation on the evolution of 
massive stars, for which predictions are now available from Heger, 
Langer & Woosley (2000), Heger & Langer (2000) and Meynet 
& Maeder (2000). These authors found that it makes a star of a 
given mass more luminous and longer-lived. The models also pre- 
dict changes in the surface chemical composition in the sense that 
He and N become enriched whilst C and O are depleted. These 
changes are stronger for higher masses and/or initial rotational ve- 
locities. 

To our knowledge, the only study which compares rotational 
evolutionary models with empirical data on high-mass stars in close 
binaries is that of Hilditch (2004). There is a known discrepancy be- 
tween masses inferred from the location in the HR diagram relative 
to the model evolution tracks, and dynamically determined masses 
(Herrero et al. 1992). Hilditch found that the discrepancy did not 



disappear when rotating stellar models were used. We begin our 
comparison against V380Cyg with the models of Schaller et al. 
(1992) and Claret (1995) for solar metallicity (Z = 0.02), moder- 
ate overshooting (qov = 0.2) and no rotational effects. Both good 
and bad agreement is seen for the Schaller models (Fig.[9l left pan- 
els): star A is overluminous for its mass, whereas the predictions for 
star B perfectly match its position. From the isochrones (green dot- 
ted lines) we estimate the age of the stars to be r = 9.5 ± 0.5 Myr. 
A similar situation occurs for the Claret models (Fig.[9l centre pan- 
els) for a metal abundance of Z = 0.02, whilst the models for 
Z = 0.01 fail to match both stars. In the TcS — logg diagram 
the situation is the same: the models can match star B but strongly 
disagree with the properties of star A. 

Fig.[9] (right panels) compares evolutionary models including 
rotational effects with the observed properties of V380Cyg. We 
interpolated in the tables published by Geneva groufjf] (Ekstrom 
et al. 2008). The models are parametrised by the ratio of the ro- 
tational velocity to the critical velocity, Jl/ficrit- The models for 
Jl/ficrit =0.1 (green tracks), and 0.5 (blue tracks) both fail to 
match the luminosity of star A. The discrepancy amounts to about 

* |http : //obswww. unige ■ ch/Recherche/evol/tables_rotzams71 
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1.5 Mq. Star B is also slightly above evolutionary tracks for its 
mass (dotted lines), but agrees within the errorbars. The Tea — log g 
diagram reveals the same story. It is interesting to note that in ro- 
tational evolution models the position of star A is around a critical 
evolutionary stage. Since rotational mixing extends the MS lifetime 
of high-mass stars, depending on initial rotational velocity, star A is 
either in the hydrogen-shell-burning phase or just at the 'blue hook' 
at the TAMS. 

Further insight into V380 Cyg's evolution comes from its pho- 
tospheric chemical composition. We interpolated in the Ekstrom 
et al. (2008) tables to the mass of star A. For f2/r2crit = 0.10 
(Vini ~ 40 kms^^) the increase in the photospheric helium abun- 
dance at the end of core H-buming is 0.5%, whilst N/C = 0.69 and 
N/O = 0.24. Increasing to il/ficrit = 0.30 gives a helium enhance- 
ment of 6%, and a further increase in N at the expense of C and O 
(N/C = 1.12 and N/O = 0.40). Our analysis of the helium abun- 
dance in star A's photosphere (Sec. l6.3t . gives a solar value, and 
we concluded that no enrichment is observed. Also, abundances of 
CNO elements (Sec.[64) give N/C = 0.25 and N/O = 0.11 which 
is exactly the initial (ZAMS) values for the models. The photo- 
spheric chemical composition appears to be untouched even though 
the star has reached the end of the core H-burning phase. Models 
for n/f^crit = 0.10, 0.30 and 0.50 at TAMS would have rotational 
velocities of about 20, 80 and 160 kms^^, respectively. Since star 
A's rotational velocity is now ~I00kms~^ we interpolated in the 
tables to find that at the TAMS a 13.13 M© model predicts N/C 
~ 0.40 and N/O ~ 0.14, while initial values for this model are 
(N/C)i„i ~ 0.30 and (N/0)i„i ~ 0.11. The later ratio is exactly 
what we found for star A, whilst no changes are observed in N/C 
ratio. We conclude that no changes in the abundances of the CNO 
elements were detected, which is contrary to model predictions. 

Why might this be? The evolution of stars in close bina- 
ries can be modified by several effects. In particular, tides affect 
their rotational velocities and so the efficiency of rotational mix- 
ing. This has been studied theoretically using detailed calculations 
by De Mink et al. (2009), but unfortunately only for masses larger 
than V380Cyg A. De Mink et al. found large N enhancements for 
the shortest-period binaries, for which rotational velocities are the 
highest, and for high-mass primary stars which are close to fill- 
ing their Roche lobes. But, for the relatively long orbital period of 
V380Cyg, 12.4 d, surface abundance changes would be small (and 
probably too small for us to have detected). We hope that future 
models will extend to the parameter space occupied by V380 Cyg, 
allowing a direct comparison between observation and theory. 



9 SUMMARY 

We have presented a study of the eccentric eclipsing and double- 
lined spectroscopic binary V380 Cyg. The system is particularly in- 
teresting because it contains an evolved star which is at the TAMS, 
so is ideal for probing theoretical evolution models which include 
rotationally induced mixing. 

New high-resolution and high-S/N spectra were secured at 
several telescopes, and a spectral disentangling technique was ap- 
plied to obtain the individual spectra of the components, and to 
(simultaneously) derive the spectroscopic orbit. Although the sec- 
ondary star contributes only a very small fraction of the system 
light (6%), we were able to attain a much higher precision than 
previously achieved. Coupled with a reanalysis of the existing light 
curves of this system, we have determined the masses and radii 
of the two stars to accuracies of 1-2%. The light ratios found in 



the light curve analysis also allowed the disentangled component 
spectra to be accurately renormalised to the continuum level. These 
spectra were used to derive the effective temperatures of the two 
stars using two independent methods, helped by the accurate sur- 
face gravity values obtained from their known masses and radii. 

The spectrum of V380 Cyg is dominated by the primary star, 
and its disentangled spectrum has a S/N approaching 1000 per pixel 
and a wavelength coverage of 3900-7000 A. We performed a de- 
tailed abundance analysis by fitting non-LTE spectra to the ob- 
served line profiles. We found no abundance excesses for helium 
or metallic species, and in fact the abundances are those of a typ- 
ical B star. The helium abundance analysis showed a broad range 
of values from different lines, which we attribute to the failure of 
calculations for the lines formed higher in the photosphere due to 
the asphericity of the star. This problem deserves to be investigated 
further It also remains possible that helium stratification might be 
the cause. 

We have compared the bulk properties of V380Cyg to theo- 
retical calculations from stellar models both with and without rota- 
tional effects. No set of models fully satisfies the observations, as 
the primary is always overluminous for its measured mass. We are 
still faced with well-known discrepancy between masses inferred 
from evolutionary models and those which have been dynamically 
measured. The mass discrepancy amounts to about 1.5 M0 for 
models including rotational evolution, and about 2.5 M0 for mod- 
els which include only core overshooting, even though we have 
revised the mass of V380 Cyg A upwards by a substantial amount. 
The surface chemical composition of high-mass stars should be an 
important probe of theoretical models. Recent model calculations 
which account for rotationally induced mixing predict changes in 
the abundances of elements involved in the CNO cycle in high- 
mass stars: helium and nitrogen should be enriched and carbon and 
oxygen depleted (strongly depending on the mass and initial rota- 
tional velocity). For a 13 Mq star all available models predicted 
changes at the end of core hydrogen burning, but the abundance 
pattern we derived in this work resembles that of a star at the 
ZAMS. We can conclude only that a detailed analysis performed 
for V380 Cyg A does not corroborate theoretical predictions of the 
chemical evolution of high-mass stars. We are currently engaged in 
extending the sample of close binaries for which such comparisons 
can be made. 
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APPENDIX A: OBSERVING LOGS FOR THE DATA 
PRESENTED IN THIS WORK 

The tables in this section contain observing logs for each of the 
fours sets of spectroscopic observations used in this work. In each 
case the orbital phases have been calculated with ephemeris derived 
by Guinan et al. (2000). 
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Table Al. Observing log for the Ondfejov red spectra of V380 Cyg 



Table A2. Observing log for the Ondfejov blue spectra of V380 Cyg. 
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Table A4. Observing log of the NOT and CAHA spectra of V380 Cyg. 
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